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Abstract 

The Schrodinger equation is solved for an A-nucleon system using an expansion of the wave 
function in nonsymmetrized hyperspherical harmonics. Our approach is both an extension and a 
modification of the formalism developed by Gattobigio et al. The extension consists in the 

inclusion of spin and isospin degrees of freedom such that a calculation with more realistic NN 
potential models becomes possible, whereas the modification allows a much simpler determination 
of the fermionic ground state. The approach is applied to four- and six-body nuclei (^He, ^Li) with 
various NN potential models. It is shown that the results for ground-state energy and radius agree 
well with those from the literature. 



1 



I. INTRODUCTION 



Expansions of the wave function of an A-nucleon system on a specific basis set is a 
common tool in few-nucleon ab initio calculations (for an overview see Ref. jsl). One of the 
challenges in such calculations is the implementation of the proper permutation symmetry, 
i.e. of the antisymmetry for a fermion system like the atomic nucleus. A frequently used 
basis are the hyperspherical harmonics (HH). For a growing particle number an effective 
HH symmetrization method had been developed in Ref. {4]. The method, however, requires 
considerable computational resources with respect to memory and cpu time. Recently a 
different HH approach has been proposed in order to construct wave functions with a proper 



permutation symmetry 



^\. It consists in the use of a nonsymmetrized HH (NSHH) basis 
thus avoiding the HH symmetrization procedure. Since the Hamiltonian commutes with the 
permutation operator, all non degenerate eigenstates of the Hamiltonian have a well defined 
permutation symmetry. The disadvantage is that additional effort has to be devoted to the 
determination of the respective symmetry. On the other hand, in view of an application of 
the HH expansion method beyond A = 4, the use of NSHH might be advantageous. 

In this work we extend the NSHH to consider spin and isospin degrees of freedom, which 
allows us to work with modern realistic NN potentials. In addition we show how the an- 
tisymmetric ground state can be identified in a very simple and effective way. We apply 
our technique to calculate ground-state energies and radii of ^He and ^Li with various NN 
potential models. In case of ^He we also use the modern realistic AV18 potential jsj. 

Our paper is organized as follows. The HH basis with inclusion of spin and isospin 
degrees of freedom is discussed in Section [TTl The formalism of Ref. for an efficient 

use of a nonsymmetrized HH basis is laid out in the first part of Section llllt whereas our 
modifications of the method are outlined in the second part of Section lllli Our results for 
^He and ®Li are discussed in Section llVt which also includes a brief summary. 
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II. THE HH BASIS 



To describe the HH basis, we first introduce the Jacobi coordinates. Here we use them 
in reversed order and for particles with equal mass, 



A-% 
A + 1 -i 



z = l,...,iV, 



(1) 



where A is the number of particles, = A — 1, and is the position vector for the i-th 
particle. To each Jacobi vector 77^ an angular momentum operator Zj is associated. For the 
total orbital angular momentum operator up to the ri-th coordinate (n = 2, . . . , A^) one has 

-C'n = -E'n-l + ■ (2) 

with Li = Zi. The hyperradial and the hyperangular variables pn and 9n are defined as 



sm 6'n = — , Pn 

Pn 



J2vh n = 2,...,N. 
\ i=i 



(3) 



The 3A^ variables {r}-^^, . . . , ry^} are thus replaced by {p, 0,]^}, i.e. one hyperradial coordinate 
p = Pj^ and 3A^ — 1 hyperangular coordinates fijv = {171, . . . , t/tv, 62, ... , 0^}. The following 
relation, 



\ 



1 ^ 



(4) 



j>i=i 



shows that the hyperradius is symmetric under permutation of particles. 



Using the above set of coordinates, the Laplace operator for n =1,. . . , can be rewritten 



as 



1 •^.3n-l^_1^2 



i=l 



pin-1 Qp Qp p2 

92 M-l d 



dp^ 



—^2 

p dp p2 



(5) 



where K"^ is a generalization of the usual angular momentum operator = l\ and is 
called the grand angular momentum operator. Its explicit expression for n > 2 is given by 



d'^ 3n - 6 - (3n - 2) cos(26'„) d 



del 



sin(2^^ 



cos2 6, 



sin2 9n " 



(6) 
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The hyperspherical harmonic functions y^K] {^n) are the eigenf unctions of the grand 
angular momentum operator: 



(7) 



The HH can be expressed in terms of the spherical harmonics Yim.{fl) and of the Jacobi 



polynomials P^' {z): 



N 



^ {limil2m2\L2M2){L2M2hm3\L3Ms) . . . 

m.i,...,mjv 

AT 

: {LN-iMN-ilNmN\LNMN) Y\ iVj] 



N 



J=2 

^pfe+l/2],[i^,_, + (3i-5)/2] (cos (2^,-))" 

The coefficients Mj{Kj \ LjKj^i) are normalization coefficients given by 



(8) 



j-i) 



l^j\{2Kj + 3j - 2)r fij + Kj^i + I 



-, 1/2 



r(/ij + ij + |)r + Kj_^ 



(3j-2) 



(9) 



where the numbers Hj are non-negative integers, and Kj = Kj-i + + Ij. 

The basis functions consist not only of the hyperangular part, but also of the hyperradial 
part and of the spin and isospin part, 

m = l-JZr^yiKu) ® \^^n ' 

where i enumerates the basis state $j, and TZ^ are the hyperradial functions numbered by 
the index r (we use Laguerre polynomials). The hyperspherical harmonic functions y[K] are 
identified by the following set of quantum numbers 



[K] = {K^,K^.u...,K2,Ln,L 



N-1, 



L2, In, ^N-i, ■ ■ ■ ihi h} 



(10) 



and X[5] ) X[T\ ^Pi^ ^^"^ isospin basis states, respectively, defined by the sets of 

quantum numbers 



\S] — {S S A-1-, ■ ■ ■ ■, S2} 

[T] = {Ta,Ta-i,...,T2;Ta,.}. 



(11) 
(12) 



The spin quantum number Sj is the value of the coupled spin momenta of particles 1 to 
j. The isospin states are characterized in a similar way with the addition of Ta,z, the third 
component of the total isospin Ta- 

In order to consider potentials that depend on spin and isospin, one has to take as basis 
states the complete functions and so one has to consider the whole set of quantum 

numbers {[K], [S], [T]} limited by a maximal value Kmax of Kn- Two possibilities arise: 
if the potential is central, then not only the total angular momentum J, the total isospin 
T = Ta (isospin mixing neglected) and its projection Ta,z are good quantum numbers for 
the eigenstates of H, but also the total orbital angular momentum L = L^r and therefore the 
total spin S = Sa (in fact J becomes obsolete). In this case the total number of basis states 
is given by the number of radial functions times the number of hyperspherical harmonic 
functions times the number of spin-isospin functions, subject to the condition that S, L, T 
and Ta.z have the desired values. On the other hand, if the potential is noncentral, like for 
example the realistic potential AV18, L and S are no longer good quantum numbers, leaving 
only J, T and Ta^z as constraints in the construction of the basis, and of course leading to a 
higher number of basis states. The size of the basis is the limiting factor in the applicability 
of the method. Examples for the number of basis states are presented in Section IIVI 



III. EXPANSIONS WITH NONSYMMETRIZED HH FUNCTIONS 

' n order to describe the nonsymmetrized HH (NSHH) method introduced by Gattobigio et 
al. l|,|2j, we first discuss a few aspects concerning permutation symmetry. The HH functions 
defined above do not have well-defined permutational symmetries. A particle permutation 
changes the definition of the Jacobi coordinates as well as the couplings between the different 
angular momenta. Consequently, the effect of particle permutation on the HH functions is 
rather complicated, and the matrices representing the permutation operators Pij on the 
HH basis are arbitrary matrices (not diagonal or block-diagonal matrices), and have to be 
calculated numerically. In order to reduce the evaluation of the matrix elements of the 
potential acting between particles i and j into a one dimensional integral one has to use 
Jacobi vectors such that r]j^_i = \J\{ri — rj). This transformation can be realized through 
the permutations A./i-i and Pj^A acting on the HH basis states. The matrix representing 
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the permutation operator Pjj can be written as 

^Ik,S,T][K',S',T'] = {^[K,S,T]{^n)\^[K' ,S' ,t']{^n)) , (13) 

where the set of quantum numbers [K, S,T] = { [K] , [S] , [T] } identify the orbital, spin and 
isospin quantum numbers of the basis functions, fijv is the set of all hyperangular variables 
defined above and fi^ indicates the set of hyperangular variables with the particles i and 
j interchanged. All permutations can be expressed as a product of a certain number of 
transpositions (depending on the particles indexes i and j) in which only adjacent particles 
are interchanged, 

4- = n = 4- • • K ' (14) 

1=1 

where the operator exchanges particles k and k + 1. This can be interpreted as follows: 
each particle carries a number and occupies originally a a box with the same number. To 
bring particle number i in the box j and vice versa, one can act with a certain number Uij 
of permutations that exchange only particles in adjacent boxes. Thanks to the properties of 
the HHs, the matrices 

^[k^!t][k',s',t'] = {^Ik,s,t]{^''''^^^)\^[k',s',t']{^n)) (15) 

representing the operators Pk are block diagonal and can be easily calculated. This provides 
an easy and fast way to calculate the product of the B matrices on state vectors. 

We remark that the use of the NSHH has been developed in order to avoid the need for the 
symmetrization of the basis functions, which, as already pointed out in the introduction, 
requires non negligible computational resources. Though the number of basis functions 
for equal values of iTmax is in general considerably larger in the NSHH method than with 
symmetrized functions, it might be still advantageous to use nonsymmetrized functions. 



In Ref. |2|] bound states have been calculated using the Volkov potential jo], which is 
central and independent on spin and isospin. In this case the spatial wave function can be 
treated separately from the spin-isospin part. The Hamiltonian is represented on the HH and 
diagonalized. Then the symmetry of the eigenstates is analyzed and the spin-isospin part 
with the correct symmetry is multiplied in order to obtain an antisymmetric wave function. 
It is clear that following this procedure the number of HH basis states, and therefore the 
size of the Hamiltonian matrix to be diagonalized, is relatively low. However, this procedure 
only applies to potentials that do not depend on spin and isospin. 
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In the following we describe how to modify the method outlined above in order to be able 
to treat also realistic potentials. In this case we need to deal with basis functions in the J- 
coupling as described in Section HIl This leads to an increase in the size of the Hamiltonian 
matrix and therefore in the number of eigenfunctions that need to be analyzed in terms 
of symmetry. In order to speed up the search for the antisymmetric ground state we have 
implemented a method based on the use of the transposition class sum operator, the Casimir 
operator, of the permutation group. This method is analogous to the Lawson method Q] 
for the removal of the spurious center of mass motion in Shell Model calculations. 

The Casimir operator of the permutation group C{A) = '^f^^^i Pij and the Hamiltonian 
H commute; hence they can be diagonalized simultaneously. In general, the eigenvalues 
of the operator C{A) are not sufficient to identify the irreducible representations of the 
permutation group. However, the completely symmetric and antisymmetric representations 
correspond to the extreme eigenvalues of C{A) and are well separated from the rest of 
the spectrum. Antisymmetric states correspond to the lowest eigenvalue of C{A), while 
symmetric states correspond to the highest one, in detail one has 

C{A)<ifM = Xm^m , 

CiA)^A = -^^^Y^^A = Xa'^a , (16) 

where ^^5, "^m, ^a, correspond to symmetric, mixed-symmetry, and antisymmetric states, 
respectively, and < Xm < A5. We diagonalize the matrix 

H' = H + ^C{A) , (17) 

where 7 is a real parameter. The eigenvalues of H' are given by 

K,r = EnT + 7Ar , (18) 

where En,r {'n = 0, 1, 2, . . . , A'max(r)) are the eigenvalues of H for the symmetry F 
(r = S,M,A). We denote by S the lowest eigenvalue of the Hamiltonian H {£ = 
min{i?o,5, Eq,Mi -E'o.a}), which can be found by performing the calculation with 7 = 0. 

To calculate the lowest antisymmetric eigenvalue of i.e. -Eq.a, we choose 7 > large 
enough so that E'q ^ is by far the lowest eigenvalue of H' . Thus one imposes the relations 

E'^^A = Eo,A + lXA<En,T + l\T Vn, Vr = ,S,M, (19) 
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which corresponds to 

7 > ' f . (20) 

Assuming that Eq a < 0, and using the fact that Ar — Ayi > A, the following condition is 
sufficient 

7>^. (21) 

Thus, with a proper value of 7 the lowest eigenstate of the Hamiltonian matrix H' is the 
physical antisymmetric wave-function, and the correct value of the ground-state energy Eq^a 
is obtained by subtracting 7AA from E'^^^. Of course with a proper choice of 7 one can 
also calculate excited states. To this end, there is no need to calculate more than a few 
eigenstates of the Hamiltonian matrix, and one can use the Lanczos algorithm. 

IV. DISCUSSION OF RESULTS 

Within the present formalism we have calculated ground-state energies and radii of 
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He and ®Li with various NN potentials using the following models: Volkov jo] (central), 
MTI/III 8|] and Minnesota (MN) 9|] (central spin-isospin dependent); in addition for '^He 
we have used AV4' 10] (central spin-isospin dependent) and the realistic AV18 jsj. In all 
our calculations the Coulomb force is included and the isospin mixing is neglected. Since for 
Volkov, MN, and MTI/HI interactions there is no unique parameter setting we define these 
potential models in Table [B 

To illustrate the typical size of the HH basis, in Table HTl we list its dimension for increasing 
-^max- It is evident that the number of basis functions grows very rapidly. Comparing for 
example the dimensions with i^max = 10, one sees that for ^He the number of HH states 
grows from about 1000 for a central interaction to about 5000 for a noncentral potential. For 
^Li one already starts with about 500000 states with a central force which is then increased 
to about 7 million states with a noncentral potential. Including the hyperradial part the 
size of the Hamiltonian matrix increases further. In fact the numbers in the table have to 
be multiplied by a factor 20 which is the number of Laguerre polynomials that we consider 
for each K. 

To improve the convergence we have used the EIHH (Effective Interaction Hyperspher- 



ical Harmonics) formalism 



111 ]. To distinguish between the calculations without and with 



effective interaction we use the acronyms NSHH and EI-NSHH, respectively. The use of 



TABLE I: Definition and parameter setting for Volkov, MN, and MTI/III NN potentials. Volkov: 
V{r) = l/iexp[-(r/;Ui)2] + V2exp[-(r//i2)^]; MN: V{r) = exp(-^ir2) + F2 exp(-/i2r2) and 
V{r) = Vi exp(— /iir^) + V3exp(— ^3r^) for NN channels with spin 5 = 0, isospin T = 1 and 
S = 1,T = 0, respectively; MTI/III: V{r) = exp(-/iir)/r + I/2 exp(-/i2r)/r {S = 0, T = 1), 
V{r) = V\ exp(— /iir)/r + V3 exp(— /i3r)/r (5" = 1, T = 0). The parameters Vi are in units of MeV 
(Volkov, MN) and of MeV-fm (MTI/III); the parameters iii are in units of fm, fm~^, and fm^"*^ for 
Volkov, MN, and MTI/III potential, respectively. 





Volkov 


MN 


MTI/III 


i 


Vi l^i 


Vi fj'i 




1 


144.86 


0.82 


200.0 


1.487 


1458.047 


3.11 


2 


-83.34 


1.60 


-91.85 


0.465 


-520.872 


1.555 


3 






-178.0 


0.639 


-635.306 


1.555 



effective interaction accelerates the convergence considerably, as can be seen in Table IIIIl 
where we list "^He ground-state energies and root mean square (RMS) radii obtained with 
the Volkov and the MTI/III potentials. One notes that the use of an effective interaction 
is not very important for the Volkov potential, whereas it has a great effect in case of the 
MTI/III interaction. Considering for example the low value of -f^max = 4 one obtains already 
results close to the converged ones for all cases, except for the MTI/III without effective 
interaction. The difference is explained by the fact that the Volkov force has a very soft core, 
whereas the MTI/III potential has a rather strong short-range repulsion, which requires HH 
functions with a rather large K]\f. Thus, in order to accelerate the convergence, we use the 
EIHH method for all the results discussed further below. In the table we list in addition 
a selection of results from the literature (note that for the MTI/III potential parameter 
settings different from ours are also in use). Comparing with our NSHH results one observes 
a very good agreement. 

In Table [IV] we list results for ^He obtained with MN and AV4' potentials. Also in this 
case we find a very good convergence and a good agreement with results from other authors. 
As already shown in Table HTl the number of HH states increases significantly using a realistic 
NN force instead of a central potential model. Nonetheless, as Table |V] illustrates, we are 
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TABLE II: Number of NSHH basis functions (including spin-isospin functions, but without hyper- 
radial part) for increasing values of -f^maxi for ^He and ^Li with central and noncentral potentials. 





^He 


6Li 


-^max 


central noncentral 


central 


noncentral 


2 


24 


54 


675 


2750 


4 


84 


264 


5400 


40025 


6 


224 


852 


30600 


315675 


8 


504 


2172 


137025 


1728950 


10 


1008 


4746 


514215 


7392960 


12 


1848 


9269 


1678950 


26377350 


14 


3168 


16776 






16 


5148 


28404 






18 


8008 


45694 






20 


12012 


79488 






22 


17472 


104988 






24 


24752 


151788 







also able to reach convergent results for ^He with the modern realistic AV18 potential. In 
the table we also make a detailed comparison with an EIHH calculation with symmetrized 
HH states. It is evident that both calculations lead essentially to the same results. Further 
results from other groups displayed in Table IVl also confirm the high precision of present-day 
few-nucleon calculations. 

Now we turn to the ^Li nucleus. In Table |VT] we present results for central potential 
models only. The reason is the non-parallel character of our present code, which prevents us 
from using a very large basis (in future we plan to work with a parallel code which should 
allow to obtain converged results for A = 6-8 with more realistic nuclear force models). The 
Table shows that, similarly as for ^He case, the results for the Volkov potential converge 
already for rather low /Tmax values. On the contrary for the MN interaction and even more 
for the MTI/III potential it would be desirable to make calculations for even higher iCmax 
values in order to obtain a more converged result for the radius. However, our ^Li results are 



10 



TABLE III: '^He ground-state energy Eq and radius trms in units of MeV and fm, respectively, 
with Volkov and MTI/III potentials with (EI-NSHH) and without (NSHH) effective interaction. 



Volkov 



MTI/III 



NSHH 



EI-NSHH 



NSHH 



EI-NSHH 



K 


Eq 


rRMS 


Eq 


'^RMS 








Eq 


rRMS 


2 


-28.579 


1.4877 


_ 


_ 


— 




_ 


— 


— 


4 


-29.281 


1.4920 


-30.414 


1.4907 


-9.679 


1 


8667 


-30.782 


1.4236 


6 


-29.811 


1.4867 


-30.417 


1.4905 


-13.744 


1 


6401 


-30.599 


1.4263 


8 


-30.160 


1.4847 


-30.449 


1.4901 


-19.923 


1 


4612 


-30.922 


1.4193 


10 


-30.276 


1.4854 


-30.407 


1.4899 


-22.807 


1 


4196 


-30.584 


1.4235 


12 


-30.363 


1.4871 


-30.422 


1.4902 


-25.808 


1 


3984 


-30.763 


1.4221 


14 


-30.390 


1.4881 


-30.416 


1.4900 


-27.253 


1 


3938 


-30.679 


1.4223 


16 


-30.405 


1.4889 


-30.417 


1.4900 


-28.385 


1 


3946 


-30.700 


1.4221 


18 










-29.061 


1 


3978 


-30.687 


1.4221 


20 










-29.558 


1 


4015 


-30.696 


1.4220 


22 










-29.856 


1 


4048 


-30.687 


1.4219 


24 










-30.092 


1 


4081 


-30.693 


1.4220 



NSHH \2\ 
HH [12,13] 
SVM [14] 
EIHH [nj 
CHH [15j 



-30.418 

-30.420 1.490 
-30.42 



-30.71 1.4222 
-30.69 1.421 



of the same level of precision as those from the literature and the agreement is satisfying. 

We summarize our work as follows. The Schrodinger equation is solved for an A-body 
nucleus using an expansion of the wave function in nonsymmetrized hyperspherical harmon- 
ics. Our approach is both an extension and a modification of the formalism developed by 



Gattobigio et al. 



21]. The extension consists in the inclusion of spin and isospin degrees of 



freedom such that a calculation with more realistic NN potential models becomes possible. 
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TABLE IV: EI-NSHH results for Eq and trms (units as in Table ^ of ^He with MN and AV4' 
potentials. 





MN 


AV4 




K 

-''-max 


Eq 


^'RMS 


Eo 




2 


-29.723 


1.4049 


-32.258 


1.3698 


4 


-30.065 


1.4139 


-32.227 


1.3840 


6 


-29.950 


1.4112 


-31.781 


1.3957 


8 


-29.981 


1.4108 


-32.579 


1.3798 


10 


-29.937 


1.4104 


-31.858 


1.3883 


12 


-29.951 


1.4107 


-32.201 


1.3866 


14 


-29.945 


1.4104 


-32.047 


1.3865 


16 


-29.946 


1.4104 


-32.068 


1.3865 


18 


-29.945 


1.4104 


-32.051 


1.3865 


20 


-29.945 


1.4104 


-32.060 


1.3865 


22 


-29.945 


1.4104 


-32.049 


1.3862 


24 


-29.945 


1.4103 


-32.054 


1.3863 


EIHH [nj 


-29.96 


1.4106 






HH [12,13] 


-29.947 


1.4105 






SVM [14] 


-29.937 








GFMC [10] 






-32.11(2) 





Different from Refs. [l|, |2| the cumbersome application of projection operators in order to 
determine the permutation symmetry of the various eigenf unctions of the Hamiltonian H is 
avoided. Instead, a potential term proportional to the Casimir operator C{A) is added, such 
that, after diagonalizing the new Hamiltonian, the state with the lowest energy is manifestly 
antisymmetric. Since one knows the eigenvalue of the additional potential term, and since H 
and C commute the ground-state energy and wave function of the original A-hodj problem 
can be obtained in a straightforward manner. In this work the approach has been used to 
calculate binding energy and radii of four- and six-body nuclei (^He, ^Li) with various NN 
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TABLE V: EI-NSHH results for Eq and trms (units as in Table |ITll) of ''He with AV18 potential. 





present work 


Reference [16] 


K 

-''-max 


Eq 








2 


-24.640 


1.5063 


_ 


— 


4 


-26.124 


1.5111 


— 


— 


6 


-25.311 


1.5061 


-25.312 


1.506 


8 


-24.999 


1.5089 


-25.000 


1.509 


10 


-24.442 


1.5197 


-24.443 


1.520 


12 


-24.491 


1.5176 


-24.492 


1.518 


14 


-24.348 


1.5184 


-24.350 


1.518 


16 


-24.313 


1.5181 


-24.315 


1.518 


18 


-24.271 


1.5177 


-24.273 


1.518 


20 


-24.266 


1.5176 


-24.268 


1.518 


22 


-24.246 


1.5170 






HH [12,i3j 


-24.22 


1.512 






FY [17] 


-24.23 








FY [18] 


-24.22 


1.516 






AGS [19] 


-24.24 









potential models and it has been shown that the results agree with those from the literature. 
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